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Abstract 

We characterize the tail asymptotics for a pure diffusion process (Theo- 
rer dl.21 . and a time-dependent local volatility model (Theoren il.311 . using 
results by Doss |D77| and NorrisfcStroock [NS91| respectively. The asymp- 
totics obtained are expressed in terms of a certain geodesic distance, and in 
the general case in terms of an Energy functional. We compare these two 
quantities qualitatively and numerically, by solving the Euler-Lagrange 
equation. We also describe the large-strike implied volatility smile asymp- 
totics, using the right-tail- wing formula of BenaimfcFriz [F3F06 II . We de- 
rive a similar result for the CEV process (Theorem OJ, and a CEV pro- 
cess evaluated at independent stochastic time r(u. t '(('Theorem 12. lft . The 
latter is applicable to the CEV-Heston model introduced by AtlanfcLeblanc |AL05f . 
Finally, we show that if we wish to use an extended version of the Carr- 
Lee |CL04] methodology to infer the characteristic function of r(t) from an 
observed single-maturity smile under the time-changed CEV model, then 
the tails of the distribution function of r(t) must have sub-exponential 
behaviour. 
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Introduction and literature review 

Davies |Dav87| introduced a new technique, using a log Sobolev inequality, for 
obtaining sharp bounds for the heat kernel, in terms of the geodesic distance 
under the Riemmanian metric induced by the inverse of the second order co- 
efficient of the diffusion (sec Appendix). However, the operator governing the 
diffusion process had to be in the so-called divergence form 

L = V-[o(a;)V] (0.1) 

and the coefficients had to be time-independent. NorrisfcStroock |NS9 1| suc- 
ceeded in removing these restrictions, but the bounds thus obtained were now 
expressed in terms of a more general Energy functional E(t, x; u 1 y), which is the 
solution to a variational problem which is more difficult to compute than the 
Riemmanian distance function. We apply their result to a time-inhomogneous 
local volatility model in Theorcir lOl 

In the time- independent case, our approach (in one dimension) is simpler, 
and merely consists of bounding (or more precisely, "sandwiching" ) the Ito func- 
tional / : C[0,T] — > C[0,T], which maps the Brownian control process to the 
response i.e. the unique strong solution of the SDE, based on the procedure 
outlined in Doss |D77] (see section 1). We use this to characterize the tail be- 
haviour of the distribution function of the marginal distributions for a general 
onc-dimcnsional SDE, subject to certain regularity conditions, in terms of the 
geodesic distance (Thcorcir ll.21) . In general, the Ito map is only continuous 
(in the topology of uniform convergence) in the one-dimensional case (see e.g. 
chapter 2 of LyonsfcQia n|LQ02|), S ussman |Suss7 8 . Lamperti |Lam64| and page 
188 in DcmbofcZeitouni [Dem^ei93] 'l. 

BenaimfeFriz.B F06 I| used regular variation theory to show how the tail 
asymptotics of the log Stock price translate directly to the large-strike behaviour 
of the implied volatility smile. This tells us when the lim sup in Lee's mo- 
ment formula [Lee04j becomes a genuine limit. They also describe the smile 
asymptotics when all moments of the terminal Stock price exist, where the im- 
plied variance exhibits sub-linear behaviour in the wings. In their sequel paper, 
Bcnaim&Friz BF06 II develop criteria for establishing when the aforementioned 
lim sup is a limit, in the case when the moment generating function of the log 
Stock price is known. This is accomplished using Tauberian theorems, which 
look closer at the limiting behaviour of the log of the Stock price mgf around the 
critical value where the moment explosion occurs. In particular, they are able 
to characterize the behaviour of the transition densities (on logarithmic scale) 
of the Stock price and the integrated variance under the well known Heston 
model. 

We apply the right-tail-wing formula to characterize the large-strike implied 
volatility smile asymptotics for a general time-inhomogenous local volatility 
model f Theorem II .5|) . and for a CEV process H~7l For a CEV process subordi- 
nated to an independent stochastic clock, we establish a one-to-one correspon- 
dence between certain exponential moment explosions for the law of the time- 
change, and the law of the terminal Stock price. This result can be specialized to 
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the CEV-Heston stochastic volatility model discussed in AtlanfcLeblanc |AL05| . 
and for this reason we advocate using this model over the popular SABR model 
introduced by Hagan ct al. HKLW02 , because less is known about the tail be- 
haviour of the transition pdf for this parametrization (see Remark 12. 4[1 



In section 3, we discuss calibration issues for a time-changed diffusion model. 
We use Ea l2.6l to extend the Carr-Lee |CL04j methodology for reverse engineering 
the law of the time-change from the law of the composite Stock price process 
at a single time. We find that the distribution function of the time-change 
has to have sub-exponential tail behaviour or else we have to resort to analytic 
continuation. 



1 Tail asymptotics for the transition densities of 
one-dimensional SDEs 

1.1 Doss's theorem 

We recall the following result from Doss |D77] , which is given as Proposition 2.21 
on page 295 in KaratzasfeShrevelKSSI 

Theorem 1.1. (Doss(1977)) Suppose that on a certain probability space (O, T , P) 
equipped with a filtration {Tt\ which satisfies the usual conditions, we have a 
standard, one- dimensional Brownian motion (Wt,J~t,Q < t < oo). Suppose 
that a is of class C 2 (M), with bounded first and second derivative, and that b is 
Lipschitz- continuous. Then the one- dimensional stochastic differential equation 

X t =X + f {b(X s )ds+^a(X 3 )a'{X s )}ds + f a(X s )dW s (1.1) 
has a unique, strong solution; this can be written in the form 

X t (u) = u(W t (u),Y t (u)) (1.2) 

for a suitable, continuous function u : M 2 — > IR, and a process Y which solves 
an ordinary differential equation, for every u) £ O. 

Remark 1.1. By integrating Eq 2.38 on page 296 in KaratzasfcShreve |KS91| . 
we see that 

x = d(y,u) = [ —— (1.3) 
Jy cr(C) 

or u = u(x,y) = d _1 (y,x), where d(., .) is the geodesic distance function under 
the Riemmanian metric induced by the inverse of <t(x), and d^ 1 is the inverse 
of d for a fixed value of y. Y t (ui) is the solution to 

jY t {u) = f{W t {u),Y t {u)) (1.4) 
i.e. a solution to a different ODE for each value of to. 
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Theorem 1.2. Assume that b(x) + ^a'(x)a(x) = ~^a 2 (x), so that 

dX t = -^a 2 (X t )dt + a(X t )dW t (1.5) 

is the log of a driftless Stock price process. We impose the same conditions as 
in Theorem \l.l\ on b and a. We further assume that 

C\ < -5—7 — — r < C 2 

au{x,y) 

dy 

< (7i < cr(x) < a 2 < 00 (1.6) 

for all x, y, for some real constants Ci, C 2 and 01, 02. TTien we Ziawe £/ie follow- 
ing tail behaviour 

-log¥(X t >x)~ d2{ 2 ,X) (xZ-oo) (1.7) 
/or i/ie distribution function of X t . 

Proof. From the bounds in Eg 11.61 we can strengthen the inequality on p(x,y) 
on page 29 in |KS91j to 

4" < P{x,v) < 7^ (1.8) 
Then, using Eq 2.42 on page 297 in |KS91| . we see that 

Ci t + X < Y t [u) < C 2 t + X (1.9) 

Eg 1 1.31 implies that u is monotonically increasing in its second argument, so Eq 
11.21 becomes 

u(W t (u), d t + X ) < X t (u) < u(W t (u), C 2 1 + X ) (1.10) 

The bounds on a(x) in Eq 1 1.61 ensure that d maps the whole line into itself (see 
Lanxperti |Lam64| '). We can re- write this as 

cr^di + Xo.WtH) < X t {w) < d- 1 {G 2 t + X Q ,W t {uj)) (1.11) 

which implies that 

W t (w) -d{dt + X ,X ) < d(X ,X t (w)) < W t (Lu) - d(C 2 t + X ,X ) 

i.e. d(Xo,Xt) is "sandwiched" between two arithmetic Brownian motions. The 
result is then obtained by applying L'Hopital's rule to 

$ c (a;) 

x— >oo e~ ( x ~ ^) 2 /^ ^ 

where $ c (x) = 1 — $(21), and <J>(a;) is the standard Normal cdf function. From 
this we see that the drift terms 1x1/2 — d{Ci/ 2 t + Xq,Xq) are irrelevant as x 
goes large. □ 
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1.2 Tail asymptotics for a Dupire local volatility model - 
the Norris-Stroock Energy functional 

Theorem 1.3. (Stroock&Norris, 1991). Let L denote the time- dependent, sec- 
ond order differential operator 

d d d d 

L = — (a{t, x) — ) + ab(t,x)- —(ab(t,x)) +c(t,x) 

ox ox ox ox 

associated with a log Stock price process dxt = — ^o- 2 (t,xt)dt + o~(t,xt)dWt, Wt 
a standard Brownian motion, so that 

^a 2 (t 7 x) = a 
m + a(b-b) = a 

-§~x [ai) + C = ° (LU) 

where the coefficients a and b are assumed to be differentiate, and a, b and c 
are measurable functions on R x R, and a is uniformly positive and uniformly 
continuous. We further suppose that there are constants A G [1, oo) and A G 
[0, oo ) such that A -1 < a(x) < A, \b\ 2 a + \b\% + |c| < A, where \b\l = {b, ab). Let 

a G (^,1) satisfying 2 o-i ^ ^ 2 ^ e 9^ ven - Then for all A G [l,oo),A G [0,oo) 
and all T G (0,oo), we have the following asymptotic result for the transition 
density p(t, x; u; y) associated with Ea \l.iy\ 

lim inf logp{t,x;uy) + E = Q 

MSoc Q<u-t<T, x,y£,R,E(t,x;u,x,y)>AI j^ 4o-l 

logp(t,x;u,y) + E N 
limsup sup — < 



A//"oo 0<u-t<T, x,y£,R,E{t,x;u,x,y)>M logE 4a — 2 

where the Energy functional is the solution to the variational problem 



E(t,x;u,y) = \ inf / |<y(s) - a(b - fc)|„_i(s,7 s ) ds 



4 jer(t,x;u,y) J t 
2 7 er(?L,u,y) J t cr 2 (s,7(; 



j nf ,/ TITTTT^ \i (s)-a(b-b)\ 2 ds (1.15) 



where T(t,x;u,y) = {7 G C([t, u],R) : 74 = x,j u = y,and \ j s \ 2 ds < 00)}. 

Remark 1.2. W^zi/i some pam, N orris & Stroock \NS91f also proved that E is Holder 
continuous. They also mention that we have the estimate 

X\y — x\ 2 1 . , , . , X\y — x\ 2 1 , , 

^ - - A Hu-t) < E(t, X ;u,y) < + -A(«-t) (1.16) 

m terms of the coefficient bounds A and A. 
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Corollary 1.4. For fixed t,x,u, by Ea \l.lb\ we see that E(t,x]u,y) can be 
made arbitrarily large if y is sufficiently large. Thus, Theorem I implies in 
particular that for all e > 0, there exists a y such that 

e" E£lil e- E < p(t,x,u,y) < E^ +e) e - E (1.17) 
which implies the following tail behaviour on logarithmic scale 

- log p(t,x;u,y) ~ E(t,x;u,y) (\y\ / oo) (1.18) 

Note the similarity of the Energy functional to the geodesic distance dis- 
cussed in Appendix A, which we can write as 

\d 2 {x-y) = \ inf I |7(s)G-i(s,7.)<k (1-19) 

Z 4 7Sr(0,a;;l,j/) J Q 

The factor of j appears instead of the usual ^ because, in Norris and 
Stroock's definition, a(x) = ^a 2 {x) as opposed to a 2 (x). E(t,x;u,y) and 
d 2 (x,y)/2(u — t) only coincide in the time-independent case when b = b. 



Under sufficient smoothness assumptions, computing E amounts to a calcu- 
lus of variations problem, where wc have to minimize the functional 

l 

L(x,x,t)dt (1.20) 

subject to x(0) = a, x(l) = b. The optimal solution is obtained by solving the 
Euler- Lagrange equation 

d dL dL ^ 
dt dx dx 

In the time-independent case, we can write Ea 11.201 as 

L(x,y,t) = m(x)(y - fi(x)) 2 (1.22) 
where m = l/2a 2 (x), and fi(x) = a(x)(b(x) — b), so Ea II. 2 II becomes 



d_ 

dt 



m'{x){^-n{x)f -2m(x)(^-M(x)) (/(x) (1.23) 



which is a non-linear second order differential equation for x(t) (see numerics 
overleaf). 



1.3 Large-strike behaviour of the Implied volatility smile 

Theorem 11.31 is similar in flavour to Theorem 2 in Busca et al. |BBF02| . which 
shows, under certain technical conditions, that the implied volatility of Euro- 
pean put/call options have the large-strike behaviour 

x hm^ a(x, T) = (^— t f T &l {u)du\ (1.24) 



7 



Local volatility function 




-2e x 



Figure 1: The local volatility function a(x) = | — O.le 1 x + 0.4e 

Energy functional vs 1/2 d A 2(y) for local volatility model 




Figure 2: Here we have plotted the Energy functional vs ^d 2 (x,y) for the dif- 



fusion process above with x = 0,T = 1, where d(y) 
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<7(C)- 
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points is =^d (y). Remember that the higher the value, the lower the weight in 
the tail. This was computed in Maple by numerically solving the Euler- Lagrange 
equation. 
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whe 



a±(t) = lim a(t,x) (1.25) 

x — >±oo 



We now sharpen this result by taking into account how o~(x) behaves at all 
values of x, via the Energy functional. 

Theorem 1.5. For the log Stock price process associated with the operator L 
in , we have the following asymptotic behaviour for the implied volatility I{k) of 
European options of log-moneyness k, expiring at u 

I 2 (k) k 

)l-L ^ (1-26) 

k 2E(t, xq; u, xq + k) 



as k — > oo. ip : [0, oo] — > [0,2] is given by ip( x ) — 2 — 4 \\fx 2 + x — x] . This 
means that I 2 (fc) beh aves sub-linearly for large k, because the bounds in \1.16\ 
imply that the right hand side tends to zero as k — > oo. 

Proof. This follows from Eq (v) in the right-tail- wing formula in BcnaimfcFriz |BF06 I| . 

□ 



1.4 Tail asymptotics for a stopped CEV process 



In this section, we consider the constant elasticity of variance (CEV) process of 
Cox |Cox75| : 

dX t =SXt +l3 dW t (1.27) 

with j3 < 0, where Wt is a standard Brownian motion. The CEV process 
can be obtained as a simple transformation of a squared Bessel process (see 
Revuzfc Yor |R Y9 T] and Linetskv |Lin04j ). which is why it gives rise to analytical 
formulae for many derivatives of interest . Zero is an absorbing barrier for /3 € 
(— i, 0), and a reflecting barrier for /3 < — i which is why we stop the process 
at the first passage time to zero. Under this specification, bankruptcy is a 
consequence of the Stock price falling ever closer to zero, in contrast to a Poisson 
default process. A closed- form expression for the risk-neutral probability of 
absorption (i.e. bankruptcy) is given in e.g. Davydov&Linetsky[DavLinOl]. We 
refer the reader to AtlanfcLeblanc [AL05| for further details and references. 

Theorem 1.6. If X t follows the CEV process above, so that the origin specified 
as a killing boundary if (3 < —i, then the right tail of the transition density 
p(T; xq, x) for Xt has the following asymptotic behaviour 

- \ogp(T;x ,x) 



25 2 (3 2 T 



Proof. The transition density for the stopped CEV process is given in e.g. 
DavvdovfcLinetskv |DavLin01| as 

^,-2/3-3/2,1/2 -2/3 2/3 -0 -0 

Kr;*„*)= ° exp(- >J- m-frggr) (1-28) 
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where v = ^py, and I v {.) denotes the modified Bessel function of order v. 
Note that this density will integrate to less than one, because of the non-zero 
probability of absorption. Then we just use the asymptotic relation 

I v (z) ~ -±= e z (1.29) 
\/2ttz 

for z ^> \v 2 — j\, then take logs. □ 

Theorem 1.7. If X t follows the constant elasticity of variance(CEV) process 
above, then we have the following asymptotic behaviour for the implied volatility 
I(k) of standard European put/call options of log-moneyness k 

^1 „ k(x e k )-™ -5 2 f3 2 (1.30) 
k 



as k — > oo. ip : [0, oo] — * [0,2] is given by ip(x) — 2 — 4 \\J x 2 + x — x]. This 
means that I 2 (k) behaves sub-linearly, because the right hand side tends to zero 
as k — > oo. It also means that there is no term- structure of implied volatility in 
the right wing, because the right-hand side is independent of T . 

Proof. Using Eg 11.281 the density of the risk-neutral return log j^- is given by 

_ k (xoe*)-^/^ 2 _ x^ + (x e*)- 2 ? x~/ (x e^ 
m-*oe exp( 2p(j2T )l u { 52p2j , ) 

(1.31) 

Then 

(xn e k )- 213 

The result then follows from Eq (v) in the right-tail- wing formula in BenaimfcFriz p3*F06 I| . 

□ 



2 Tail asymptotics for a diffusion process evalu- 
ated at an independent stochastic clock 

2.1 The general time-changed diffusion 

We now assume that the reference entity's Stock price process St is a one- 
dimensional diffusion process (X t )t>o evaluated at a random time given by a 
increasing, right continuous process t(lj, t), which is independent of X t ; thus 

St = -^r(t) 

dX t = a(X t )dW t (2.1) 

with a(x) bounded between two positive constants, and uniformly Holder con- 
tinuous on K + . t(T) may be continuous or purely discontinuous (e.g. an increas- 
ing Levy process); see Geman.Madanfe Yor |GMY02| and ContfeTankov |CT04| 
for further discussion on this point. 
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One way to represent the distribution of the stopping time t v = inf{f > 
: Xt = y S (0,oo)} is to express its Laplace transform E x (e~ XTy l/ r <00 }) in 
terms of the increasing and decreasing solutions of the Sturm- Liouville equation 

(Af)(x)-Xf(x)=0 (2.2) 

for A > 0, (see Davydov&Linctsky DavLinOl , Borodin&Salmincn BS96 and 
KarlinfcTavlor |KT8 lj ) . Specifically, we have 

EAe- XTL hr L <oo } ) = x>L 

E x (e"^ l {TU<oo} ) = x<U (2.3) 

Under the conditions have we placed on o~(x), and 00 arc natural boundaries, 
and V'a(0+) = 0, (f>\(0+) = 00, linx,,-^ ij)\{x) = 00,1^^00 cj>\{x) = 0. 

Remark 2.1. <j)\{x) andip\(x) are called fundamental solutions of the ODE 
in Eq \2.'A They are linearly independent, and all solutions to Eq \2."A can be 
expressed as their linear combination. 

By the Feynman-Kac formula, the solution to the Cauchy problem 

dt ~ AI 

f(T,x) = ?/>a(.t) (2.4) 

has the stochastic representation 

f(t, x) = E^iMXr)) = M*) e HT ~ t} (2-5) 

and similarly for 4>\, because zero is a natural boundary. For the subordinated 
St process, by conditioning on the independent r(T) at time zero, we observe 
that 

E(e^)) =E°.So ( MSrl) =E o,s 0( MM ) (2 . 6) 

Remark 2.2. If X t is a standard Brownian motion, then the eigenfunctions are 
given by e ± ^ Xx . 



2.2 Tail asymptotics for a time-changed CEV process 



We now consider the right tail asymptotics for the transition densities of stopped 
CEV diffusion, but evaluated at an independent stochastic clock. Adding this 
stochastic volatility component means that the smile effect is less correlated to 
the probability of default, and it means that the right tail decays significantly 
slower (this is made precise below). We do not correlate the time-change with 
the CEV process, because the leverage effect is more than adequately explained 
by the power law in the CEV diffusion coefficient. 
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Theorem 2.1. For a CEV process evaluated at an independent random time 
r(t), we have the following one-to-one correspondence between exponential mo- 
ment explosions for t(u>,T) and St = X t ^ ^{oj) 

X*(T) = sup{ Ael + : Ee Ar(T) < 00} 

S\f3\ 



sup{ A G R : E( exp (-^ S^ p )) < 00} (2.7) 



if X*(T) G (0,oo), which ensures that the tails of St = -XV pr) are realistically 
fat. 

Corollary 2.2. for an?/ random variable X , a simple Chebyshev argument gives 
rise to the following result (see Benaim&Friz \BF06 Itf ) 

- lim sup lQgF(X > X) = q* = sup{g : Ee« x < 00} (2.8) 

x — >oo % 

so we have the following estimate for the tail behaviour of the distribution func- 
tion of St 

logP(^ S T P > x) 

_ lim SU p K& W L = yfW(T) (2.9) 



Remark 2.3. Corollarv \2. 2\ could be applied to extrapolate the prices of digital call 
options to large strikes, where PDE and Monte Carlo methods break down. It 
may be possible to combine this with the tail asymptotics for the integrated vari- 
ance under the Heston model reported in section 6.3 of Benaim&Friz ]BF06 IH/ . 
to prove that the lim sup here is indeed a limit. We suspect that this is the case, 
and if so, one could then invoke the right-wing-tail formula in Benaim&Friz \BF06 H/ . 
However, this appears to be difficult because we need to ignore the term in 

Ea \1.2iA in order to prove Theorem \2.1\ 

Proof, (of Theorem 12. 1[) The CEV process violates the condition on the bound- 
edness of the diffusion coefficient in section 1.1, but the first equality in Eg 12.61 
still holds. In this case, ip\(S) is given by 

V>a(S) = S*I v (yf(2\)z) (2.10) 

where 

* = *w = w\ s ~ p (2 - n) 

and v = 5pr. I v (.) denotes the modified Bessel function (see DavvdovfcLinetskv |DavLin01j ). 
Now, using the asymptotic relation in Eg 11.291 we see that for all e > 0, there 
exists a Ki(e) such that 

MS)P(S) dS < (1 + e) / e V2Xz(s) p{s) dS (2 12) 

Jk * sj2irV2\z{S) 

where p(S) is the transition density for St, and there exists a ^(e) such that 
(1 + e ) s l_ e V2A-(S) < e V2\ z (S)(i+e) for all s > R ^ e y Taking R ^ = 

y 2tt V2A z(S) 
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max(Xi(e), ^(e)) and remembering that ip\(S) is monotonically increasing in 
S, we see that 



E(Va(S t )) =oo 



E 



exp(V2Az(S , T )(l + e)) = E cxp ( 



^^(l+e)) = oo 



We proceed similarly for the lower bound. 



□ 



Example 2.1. Ifr(t) is obtained as the integral of the Cox-Ingersoll-Ross(CIR) 
square root diffusion process dv t = n(9 — v t )dt + a^/v t dW 2 , o,Vq > 0, then the 
critical value A*(T) can be computed numerically, see e.g. Proposition 3.1 in 
Andersen&Piterbara XAPlJSI or section 6.3 of Benaim&Friz \BF06 Ilf . 

Remark 2.4. For an uncorrelated SABR model where the time-change is given 
by the exponential functional of Brownian motion ' = J * e 2 ' Sa+/iS ) ds with 
oo for all A > 0, T > (see Jourdain Jour04 ), so 12.11 only 



tells us that 



At(T) 



(2.13) 



for all A. The difficulty with the SABR model lies in characterizing the tail 
behaviour of r(t). At the present time, we have only been able to establish the 
Gaussian bounds 



exp 



x 2 (l + e) 



< 



t(T) 

»(-^ > x) < exp 



x 2 (l + e) 
8t 



(2.14) 



for all e > 0, if x is sufficiently large, using Jensen's inequality for the lower 
bound which is slightly less trivial. There is a closed-form expression for the joint 
density of (log S t ,v t ) under the SABR model with [3 = l,p = 0, originally ob- 
tained by McKean (see Hagan et al. |HLW04j . and MatsumotofcYor |MatsYor05 lj . 
|MatsYor05 TTj ). 



Remark 2.5. In the case when Ee Ar ( T ' < oo for all A > 0, it may should be 
possible to derive a sharper result that \2.1Y using e.g. Kasahara's Tauberian 
theorem (see Bingham et al. lBGT87Jj ) to relate the tail behaviour of r(T) and 

S -I3 — 

the tail behaviour of "jpp However, this case is of little practical use in Math- 
ematical finance, because we would typically add a stochastic volatility /time- 
change component to fatten the tails of the risk-neutral densities, so that some 
of the moments of the terminal Stock price explode. 



Remark 2.6. Conversely, it is also possible to derive bounds on the moment gen- 
erating function, and hence the tails, for a general time-changed diffusion when 
A*(T) £ (0,oo), using Doss's theorem on pg 295 in Karatzas&ShreveJEEMl- 
However, these are somewhat cumbersome because they involve bounding con- 
stants on the drift and diffusion coefficients, so we omit the details 
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3 Calibrating a time-changed diffusion model: 
the Carr-Lee inverse problem of extracting 
the law of the time- change from a single- maturity 
smile 

We now consider the inverse calibration problem of extracting the law of the 
time-change t(T) in Eg 12.11 given the law of St = ^ r (T), when a(x) is known. 
We would typically obtain the law of St from the observed implied volatility 
smile at maturity T. This problem was considered at length for the case when 
the X t process is geometric Brownian motion (i.e. cr(x) — 1) by CarrfcLee |CL04j . 
and later by Friz&Gatheral [FrizGath05j . 

In principle, we can extend the Carr-Lee approach to arbitrary o~(x), using 
Ea l2.6l to reverse engineer the characteristic function of t(T) from the prices of 
so-called eigenfunction contracts as E(e i9T ( T >) = E(&4§4) (for 9 g R). if the 
Expectation on the right-hand side exists. We could span the complex-valued 
eigenfunction contract paying ^'"(^j with a portfolio of standard call and put 
options, using the well known decomposition of a twice differcntiable payoff 
given in e.g. the Appendix of CarrfeMadan |CM98] . We could then compute the 
distribution function of St, using Levy's inversion theorem. For r(t) continuous, 
the most important choice for a(x) is the Carr-Lee case with o~(x) = 1, because 
then the time-change agrees with the quadratic variation of the log-returns 
process (by the Dambis-Dubins Schwarz result), which is the underlying for 
liquidly-traded volatility derivatives. Nevertheless, calibrating the law of t(T) 
in the general case is important if we wish to build a model which is consistent 
not only with the observed prices of European options, but also with barrier 
and forward starting options. The Dupire local volatility model discussed in 
section 1.4 typically falls short in this respect, because it internalizes the wrong 
dependence structure. 



When the underlying diffusion process in Eq 12.11 is a CEV process with 
a(S) = S 1+f3 , (i e (0, 1), then setting A = id for 6 e R, Eq EH becomes 



E 



Si I v {y/{2i0) z(S T )) 
S$ I v (y/(2i0) z(S )) 



E( e ^( T )) 



(3.1) 



By a very similar argument to Theorcm l2.ll we see that 



E ( e ^r(T) ) 



J S$.I v (y/{2iO)z(ST)) ' 

[Si I v {tJ[2i8)z(S )) _ 
E exp^OS^-V 3 )) 





exp( 



{i + j)Ve 

1/31 




))p(T;S ,S)dS < oo 



(3.2) 
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However, the Lebesgue integral in the last line will not exist if the integrand is 
not absolutely intcgrable, so the restriction that 



J exp (^- (S- - V))p(T; So, S)dS = E 



exp(^(5^-V)) 



< oo 



(3.3) 



provides a necessary condition on the tail behaviour of St, so we can compute 
the value of E(e WT( - T 1) for all 8 e R. This means that 



logP(r>a;) 
hm sup = oo 



(3.4) 



so the tails of the distribution function of t(T) have to decay faster than any 
exponential. This condition is violated by the Heston and SABR models, and 
may well explain why the Friz-Gatheral Moore-Penrose regularization scheme in 
|FrizOath05| performs so badly for high volatility of the instantaneous variance 
under the Heston model. Our extended Carr-Lee methodology fails in these 
cases because, for certain values of 8 in Eq 13.21 we go outside the strip of 
regularity of the characteristic function of (<Sy' 3 — where it cannot 

be represented as a Fourier integral. It is easily shown that same restriction 
applies for the original Carr-Lee framework with a(x) = 1. The only other 
way to circumvent this problem would be to compute ^f(6) = E(e ieTl - 7 " 1 ) on a 
domain G C C, and then use analytic continuation, by computing a truncated 
approximation to the power series representation for 'I'. 
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Appendix: The Varadhan density estimate, and 
the geodesic distance 

For a diffusion process with generator i ^ a ij ( x ) ~&k~ ' wnere a ij satisfies uni- 
form Holder and ellipticity conditions, the Varadhan |Var67PDE| density esti- 
mate states that 

lim -2t logp{t, x, y) = d 2 [x, y) (B. 1) 

uniformly over x and y such that d(x,y) is bounded, where 

d 



d(x,y) = inf / \fl ~S~^ qiA————)db 

V /:/(0)=x,/(l)=» J VVZ - J 



(for Revalued / € C[0, 1]) is the shortest (or geodesic) distance from x to y 
on a Ricmmanian manifold with metric gij = (a y ) , which is the inverse of 
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a = <j<j t . Under sufficient regularity conditions, we can use the Euler-Lagrange 
equation from calculus of variations theory to show that the optimal path sat- 
isfies the geodesic equation 

^ + T ^^T-° (R2) 

where the Christoffel symbol Y™p = ^g av (g v ^i3 + 9vp,n~ 9nP,v) = \a al '(g vfl ,i3 + 
9vp.fi — 9fi0,v) is the a-th component of S^fr. Geodesies parallel transport 
their own tangent vector U = i.e. the covariant derivative U@U% = 0). This 
implies that the scalar product 

is constant along geodesic paths. The geodesic length d(x, y) is independent of 
the parametrization of the curve. However, the following functional 

1 /"* df l df j 



I(x, y) = — inf / qij — — dt 

K U > 2 /:/(0)=x,/(X)=v Jo dt dt 



1,3- 



which crops up in the Wentzcll-Freidlin theory of large deviations for SDEs(see 
Varadhan [Var84] ) . does depend on the parametrization. However, for a class 
of parametrizations all representing the same curve in K d , I(x,y) is minimized 
when the curve is parametrized so as to cover equal lengths in equal times i.e. 
I(x, y) = | J Q cdt = Tjct, d (x, y) = J Q s/cdt = s/ct so that I(x, y) = ^d 2 (x, y)/t 
(see Varadhan |Var67Probj . 



